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• • . 1 INTRODUCTION 



ABSTRACT 

We present red spectra in the region ~ A7000-8300A of the eclipsing dwarf nova IP 
Peg, with simultaneous narrow-band photometry centered at 7322A. We show that by 
placing a second star on the slit we can correct for the telluric absorption bands which 
have hitherto made the TiO features from the secondary star unusable. We use these 
TiO features to carry out a radial velocity study of the secondary star, and find this 
gives an improvement in signal-to-noise of a factor two compared with using the Nal 
doublet. In contrast with previous results, we find no apparent ellipticity in the radial 
velocity curve. As a result we revise the semi- amplitude to K2 — 331.3 ± 5.8 km s , 
and thus the primary and secondary star masses to 1.O51q;J y M and 0.33^0 05^© 
respectively. Although this is the lowest mass yet derived for the secondary star, it is 
still over-massive for its observed spectral type. However, the revised mass and radius 
bring IP Peg into line with other CVs in the mass-radius-period relationships. 

By fitting the phase resolved spectra around the TiO bands to a mean spectrum, 
we attempt to isolate the lightcurve of the secondary star. The resulting lightcurve 
has marked deviations from the expected ellipsoidal shape. The largest difference is 
at phase 0.5, and can be explained as an eclipse of the secondary star by the disc, 
indicating that the disc is optically thick when viewed at high inclination angles. 

Key words: novae, cataclysmic variables - binaries: close - binaries: spectroscopic 
- binaries: eclipsing - starsdundamental parameters - stars: individual: IP Peg 



IP Peg is a dwarf nova which brightens by ~2 mag ev- 
ery 100 days or so, with each outburst lasting ~2 weeks, 
and is important because it is the brightest known eclips- 
ing dwarf nova above the period gap. The eclipses give tight 
constraints on the geometry of the system, and in combina- 
tion with other data, such as radial velocity studies, allow 
the masses of the two stars to be determined. We should 
thus be able to compare the mass, radius and density of the 
secondary star with those of main sequence stars, and learn 
something of its structure and evolution. Unfortunately the 
radial velocity curve of the secondary star in IP Peg is prob- 
lematical, being apparently elliptical (Martin et al. 1989, 
henceforth M89). The star is believed to be in a circular 
orbit, and the apparent ellipticity is thought to be due to 
irradiation of the secondary star (alt hou gh we shall present 
an alternative explanation in Section 6.4), and an uncertain 



correction must be applied before the binary parameters can 
be determined. The result is acutely embarrassing. Smith & 
Dhillon (1998) collect together all the available masses, radii 
and spectral types for CVs. IP Peg is one of the few objects 
that does not rely on the dubious method of using the accre- 
tion disc lines to determine the white dwarf radial velocity, 
and so should give reliable parameters. However, in all the 



relationships IP Peg is a persistent offender, lying well clear 
of the supposedly less reliable points in both the mass vs. 
orbital period and mass vs. spectral type diagrams. 

A further complication is presented by the lightcurve 
of the system. High-speed photometric observations show a 
light curve dominated by a bright spot and a deep eclipse. 
The eclipse is that of the bright spot and white dwarf by the 
secondary star, superposed on a more gradual disc eclipse. 
Unfortunately the bright spot ingress overlaps with that 
of the white dwarf, making the ephemeris determination 
more problematical. Wood & Crawford (1986) were able to 
use their high-speed photometric observations to derive an 
ephemeris for IP Peg based on white dwarf egress timings. 
Later observations by Wood et al. (1989) and Wolf et al. 
(1993) show that the white dwarf egress varies markedly in 
strength and duration (10s to 300s). 

In this paper we present red spectroscopy and photom- 
etry of IP Peg. The observations and their reduction are 
explained in Section |^. After that the paper divides into 
two main threads. The first is the radial velocity study and 
its results. The extraction of the velocities is described in 
Section ^ after which we discuss the summed spectra and 
new ephemeris (Sections ^ and ^). In Section ^| we discuss 
the absence of apparent eccentricity in our data, and con- 
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elude the system has changed in some way. We use our radial 
velocity semi- amplitude to derive new system parameters in 
Section The second thread, isolating the lightcurve of the 
secondary star, is explored in Section H. 



2 OBSERVATIONS 

Photometry and spectroscopy of IP Peg were obtained over 
four nights in 1995 October from the 2.5-m INT and 1.0- 
m JKT on La Palma. Only two nights of photometry were 
tak en - simultaneous with the spectroscopy (see Section 
2^). Table contains a log of the observations. 

Data from the AAVSO show that in the latter half of 
1995 IP Peg underwent two outbursts, one in mid-September 
and one in mid-December. The data presented in this paper 
were taken in the first week of October, approximately two 
weeks after the end of the mid-September outburst, and are 
thus during quiescence. 



2.1 Photometry 

Data were taken with the JKT on 1995 October 4th & 5th 
using an EEv7-chip with a scale of 0.31 arcsec per pixel. A 
narrow band (fwhm=46A) On filter with maximum trans- 
mission at a wavelength of A7322A was used on both nights. 

The images were processed by subtracting off a single 
bias image in the case of the Oct. 4th data, and a median of 
three biases for the Oct 5th data. A single sky flat taken at 
the beginning of each night was used to flatfield the im- 
ages. The photometry was then extracted using an opti- 
mal extraction method (Naylor 1998). The photometry of 
IP Peg was divided by that of a neighbouring star (the sec- 
ond star placed on the slit during the spectroscopic run — see 
below) to remove the effects of sky transparency variations 
and, thus, the photometry is relative. The phase-folded light 
curves for both nights are shown in Figs, [i] and ^| with the 
bright-spot hump at (f> ~ —0.2 and the primary eclipse at 
6 = 0.0. 



2.2 Spectroscopy 

Data were taken on 1995 October 4th, 5th, 9th and 10th 
with the INT, using a TEK4-chip with a scale of 0.33 arcsec 
per pixel, but binned by a factor two in the spatial direction. 
A two-arcsec slit was used set at PA 58.7 degrees to obtain 
the spectrum of a neighbouring star simultaneously (the star 
used to divide the IP Peg photometry by — see ab ove) . This 
allowed us to correct for slit losses (see Section 2.4). The 
R831R grating was used on the 4th and 5th, giving ~1.22 



A per pixel; the R1200R grating was used on the 9th and 
10th, giving ~0.84 A per pixel. An arc spectrum was taken 
each time the telescope's position on the sky was changed (~ 
every hour). The nights of the 4th and 5th were generally 
clear but the nights of the 9th and 10th were affected by 
cloud, especially so the data of the 9th. An exposure time 
of 600s was used on the 4th, 5th and 9th: 400s on the 10th. 

The images for each night were first bias subtracted. 
Tungsten flats taken each night were then used, after divid- 
ing by fitted low-order polynomials, to flatfield the images. 
Extraction of the spectra was performed using an optimal 
extraction method (Home 1986, Robertson 1986), and then 



0.22 
0.20 
0.18 

M 

3 0.16 

| 0.14 

| 0.12 
OS 

0.10 

0.08 

0.06 

3 10 

3 
■o o 

"8 _K 

m o 
K -10 
K -15 
5. 



_ + _ 



40 



5.45 



5.50 5.55 
BDJD-2449990.0 



5.60 



5.65 



Figure 1. The lightcurves obtained from the photometry (dots) 
and spectroscopy (bars). The lower panel shows the percent- 
age residuals between the two, calculated as 100.0 * (-R S pcc — 
^photV^phot- The point in the top left of the upper panel repre- 
sents the average size of a photometric error. Data taken on 1995 
October 4th. 
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Figure 2. As previous figure but for the data taken on 1995 
October 5th. 



wavelength calibration was done using low-order polynomial 
fits to the appropriate arc spectra. 



2.3 Removing the effects of Earth's atmosphere 

Normally, one observes a standard star (usually a near- 
featureless O star) at regular intervals throughout the night 
to divide the program star by. Here, we use the second star 
placed on the slit and observed simultaneously with IP Peg. 
This has the advantage that it avoids the assumption of uni- 
formity of atmospheric absorption with position and time. 
Fig. ^ shows that the atmospheric absorption has been well 
corrected, including the large bands near 7600A. 

The disadvantage of this technique is that the star may 
not be as featureless as one would like, since its main selec- 
tion criterion is proximity to the target star. Comparing blue 
spectra of the star with those in Jacoby, Hunter & Chris- 
tian (1984) shows it to be a late F-star. Such stars are good 
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Table 1. Journal of observations of IP Pegasi 



uate 


1NO. OI 


Wave length 


Exposure 


Phase 


Mean Seein 


^wct. yo ) 


Prames 


Region (A.) 


Time (sec) 


Covered 


(arcsec) 


Photometry: 












4th 


81 


On filter 


120 


0.51-1.82 


2.17 


5th 


100 


On filter 


120 


0.31-1.70 


1.65 


Spectroscopy: 












4th 


18 


7100-8300 


600 


0.38-1.03 


1.85 


5th 


18 


7100-8300 


600 


0.38-1.70 


1.42 


9th 


22 


7000-7800 


600 


0.56-1.72 


1.78 


10th 


39 


7000-7800 


400 


0.83-1.79 


1.98 




7500 7600 7700 



7500 7600 7700 
WAVELENGTH (A) 



7500 7600 7700 



Figure 3. This figure shows the correction for the atmospheric 'A'-band near A7600A. The left hand panel is the second star on the slit; 
the center panel is the uncorrected IP Peg spectrum; the right hand panel after correction. All the spectra have been normalised to one 
at 7500A. Note how the drop of the TiO band heads (7550-7650A) and the KI lines (7650 and 7700A) become visible after correction. 



atmospheric calibrators at low resolution (see Mason et al 
1995), but at the intermediate resolutions we used weak ab- 
sorption features are present. In Fig. |j] we have averaged all 
our higher-resolution divided spectra with no velocity shift, 
and can see these features as "emission" lines superimposed 
on a velocity smeared M-star spectrum (since we have not 
corrected to the rest frame of the secondary star). The fea- 
tures are weak (normally around a couple of percent, with 



a few at 10 percent) and as we shall show in Sections ^ and 
^j, have a negligible effect our radial velocity study and final 
summed spectrum. 



2.4 Relative spectrophotometry 

During our spectroscopic run with the INT, two nights 
of photometric data (simultaneous with the spectroscopic 
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Figure 4. The result of a simple mean of the spectra which 
result from dividing the IP Peg spectrum by the second star on 
the slit. Only the higher-resolution spectra (1995 Oct 9 and 10) 
have been used. 



data) were taken with the JKT to test the two-star spec- 
trophotometry technique. The wavelength region of the 
spectroscopy corresponding to the photometry was isolated 
and summed according to the weights of the On filter. After 
summing, the spectrophotometry of IP Peg was divided by 
the similarly summed flux from the second star on the slit. 

Figs. |l[ and ^ show the photometric data over-plotted 
on the spectroscopic data. The lower panels in the figures 
show the percentage residuals between the two, calculated as 
100.0 * (-Rspcc — i?phot)/i? P hot- Where more than one photo- 
metric point is coincident with a spectroscopic one, a simple 
average of them is taken. For the spectroscopy point covering 
primary-eclipse egress in the data of 1995 October 4th, the 
exposure was too large long sample the fast rate of change 
in brightness properly. This resulted in the relatively large 
residual seen here. Ignoring this point in the analysis, the 
average residual size for the nights of the 4th and 5th are 
3.28% and 2.80%, respectively. That is, the difference be- 
tween relative photometry and relative spectrophotometry 
is only at the few percent level and, for the most part, lies 
within the error of the photometry. Thus the real difference 
is probably well within the measurement errors. This shows 
that performing relative spectroscopy by the inclusion of a 
second star on the slit works well. 

This result is apparently in contradiction with the work 
of Webb et al. (2000). They performed a similar experi- 
ment to ours, obtaining spectroscopic observations of the 
low-mass X-ray binary J0422+32, again with a second star 
on the slit. When they compared the relative photometry 
from their spectroscopic observations with that from simul- 
taneous photometry they found differences of order 15 per- 
cent. The author in common between these two papers has 
no explanation as to why this should be so, though we note 
that IP Peg is much brighter, and that a broadband filter 
was used for J0422+32. 



3 A RADIAL VELOCITY STUDY 

We used an iterative procedure to determine the radial ve- 
locity curve of IP Peg, in which a summed spectrum (cre- 
ated after each iteration) is used as the template. This tech- 
nique is similar to that used by Mukai & Charles (1988) and 
Home, Welsh & Wade (1993). First, the mid- white dwarf 
egress timings of Wood et al. (1989) and Wolf et al. (1993) 
were fitted with a linear ephemeris, and this ephemeris was 
initially used to phase all of our spectra. The radial velocity 
semi-amplitude (K2) value from M89 was then used to as- 
sign a stellar radial velocity component to each spectrum for 
a given phase cj> by assuming a circular orbit for the red star 
of the form K2 sin(27r</>;), where (/>, is the phase of spectrum 
i. Each spectrum was then Doppler-shifted by its orbital 
radial velocity; all such shifted spectra were then co-added 
into a grand-sum spectrum. This grand-sum spectrum (less 
the spectrum of interest — see later) was then used as a 
template spectrum against which to cross-correlate all the 
individual spectra to obtain a radial velocity study. 

The velocity extraction process is run iteratively. The 
velocities from the initial set of cross-correlations are fit- 
ted with a sine wave to return a new K2 -value and cor- 
rection to phase zero; each individual spectrum is then as- 
signed a new phase and new radial velocity from this sec- 
ond generation ephemeris and iG-value; each spectrum is 
then shifted by these new velocities and combined to make 
a new orbit-averaged template spectrum; the spectra are 
then cross-correlated against this new average, new veloci- 
ties extracted and so on, and the whole process repeats until 
convergence. 

The process was considered to have converged when the 
correction to phase zero was less than 0.001 of the orbital pe- 
riod, corresponding to a correction of less than 14 seconds. 
This value was chosen because tests showed that beyond 
this, no significant improvement within the errors of the fit 
was achieved. Convergence typically occurred within ten it- 
erations. A small number of spectra returned dubious ve- 
locities; these spectra were of lower signal-to-noise (mainly 
those of 1995 October 9th, due to occasional cloud), and 
were removed from the analysis. 

Note that to avoid any danger of 'self cross-correlation, 
the spectrum being correlated is not used in creating the 
average-spectrum template against which it is being cross- 
correlated. This results in a new template being created for 
each individual spectrum, but the large number of spectra 
which adequately sample the orbit ensures that the template 
spectrum is, to all intents and purposes, identical through- 
out. 

The first two nights of data were taken at a slightly 
lower resolution than the last two, so the data were split 
into high and low resolution data sets and each one treated 
separately. For both sets of data, two regions were used in 
the cross-correlation but only one region was used at a time. 
These regions were a 115A-wide region from A7105-7220A, 
and a 170A-wide region from A7560-7730A. The slightly 
different wavelength coverage between the two data sets al- 
lowed the shorter wavelength cross-correlation region in the 
slightly higher resolution data set to be extended blue- wards 
to A7000A. 

Table |^ presents the results of this study for the two 
regions used in each data set. A weighted mean (labelled 
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Table 2. Radial velocity data using an orbital-average as template 



A 


Notes 


N 


K 2 


To-2449998.0 2 


f c 


fe 


F 


P 1 


7105- 


-7220 


TiO 


31 


353.8 ± 14.1 


0.51609(86) 


32.14 


31.67 


0.404 


0.671 


7560- 


-7730 


TiO,K 


32 


325.6 ± 14.3 


0.51722(95) 


33.92 


33.41 


0.424 


0.658 


7000- 


-7220 


TiO 


59 


310.6 ± 9.4 


0.51841(73) 


32.42 


31.78 


1.183 


0.314 


7560- 


-7730 


TiO,K 


59 


344.3 ± 9.5 


0.51612(64) 


33.01 


32.88 


0.220 


0.803 






wT, 




331.3 ± 5.8 


0.51692(39) 










7670- 


-8319 


M89 region 


32 


327.7 ± 16.1 




30.11 


30.11 


0.00 


1.00 



1 Orbit is consistent with being circular for p > 0.05. 

2 Numbers in brackets are the errors on the same number of figures immediately to their left. 



uiE in the table) of these results provides us with a new 
mid- white dwarf eclipse (not egress) timing and K%. The 
individual radial velocity curves are shown in Fig. ^ with a 
circular orbit fit using the above weighted mean values over 
plotted as a full line; the individual orbital fits of Table @ 
are over-plotted as dashed lines. The quoted error is the 68% 
confidence interval ('ler') calculated for two free parameters 
(Lampton, Margon & Bowyer 1976). To derive this error 
the radial velocity points were given equal weight, and a xt 
assumed for the best fit. 

We compared the signal-to-noise we obtained using the 
TiO bands with that from the more traditional Nal dou- 
blet. The gain using the new technique is about a factor 
two in signal-to-noise if the results from the two TiO bands 
are averaged. We also tried cross-correlating each spectrum 
against an M-star template instead of the mean IP Peg spec- 
trum. There was an increase in noise for the fits to these 
data, presumably because the mean spectrum gives a better 
match to the spectral type than any template star. 

Finally we wished to ascertain whether the weak F-star 
features introduced by our technique for removing telluric 
absorption (see Section 2.3) had affected our values for the 
radial velocity. Specifically, there was a concern that it would 
introduce a signal at zero velocity shift, which might reduce 
the semi- amplitude. During the run we took several obser- 
vations of early-type stars to check our correction technique. 
We chose a pair of observations where an O-star corrected 
the observation of the second star on our slit well, yielding a 
telluric absorption corrected spectrum of this star. We then 
divided our corrected spectra of IP Peg by this spectrum, 
thus doubling the strength of the F-star features. We then 
repeated our cross correlation analysis for Oct. 4 and 5. The 
derived semi-amplitude changed by less than the error (+7 
km s" 1 for 7105-7220A and -4 km s _1 for 7560-7730A). It 
should be noted that this is an over-estimate of how much 
the F-star affects our study, since we will also have intro- 
duced noise, and poor atmospheric correction, all of which 
might tend to decrease the semi-amplitude. Despite this, it 
is clear that the second star on the slit makes a negligible 
contribution to the cross-correlation function, for reasons we 
explain in the next Section. 



4 THE SUMMED SPECTRUM 

Fig. ^| shows the final Doppler-shifted weighted-mean spec- 
trum. The reasons for our choice of cross-correlation bands 
should now be clear. The first corresponds to the TiO band 
heads around 7100A, the second to the TiO band heads 
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Figure 6. The Doppler-shifted, averaged spectra of IP Peg taken 
during 1995 October. LR is the lower resolution data taken on 
the nights of the 4th & 5th (upper spectrum); HR is the higher 
resolution data taken on the 9th & 10th (lower spectrum). Lower 
limit symbols mark the boundaries of the cross-correlation regions 
used in the radial velocity study. 



and KI lines (A7764.9, 7699. OA) near 760uA. The final fea- 
ture obviously from the secondary star is the Nal doublet 
(A8183.3, 8194.8A). 

The emission features are harder to identify, so to aid 
this process we show the sums of spectra in fourphase ranges 
near phase zero in Fig. Q In this, and Fig. H we can see 
emission features around 7770A, 780uA, 7870A, 7930A and 
8200 A. These features are all broad, and weaken or disap- 
pear when the disc is eclipsed in the phase 0.95-0.05 spec- 
trum of Fig. [?]. This suggests a disc origin for all these lines. 
The easiest to identify is the emission around 8200A, which 
is where the Paschen series converges, but we can also expect 
disc emission from the Nal doublet, and perhaps Call 8203A 
and maybe even Hell 8237A. Friend et al. (1988) identify 
the OI triplet at AA7771.9, 7774.2 and 7775.4A, which co- 
incides with our 7770A feature. This leaves three features 
still unexplained, of which only the 7800A feature appears 
in M89. 

In principle there is a contribution to this spectrum 
from the F-star used for telluric correction. However, as 
stated in Section 



2.3 the strongest features are around 10 



percent of the continuum level. As the summed spectrum is 
velocity shifted, these features are smeared out over 600km 
s _1 , or 15 pixels, making them less than 1 percent in depth. 
The lines will be similarly smeared in the template spectra 
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Figure 5. Extracted IP Peg radial velocity curves. (1) The region 7560-7330A for Oct. 4th & 5th; (2) the region 7105-7220A for Oct. 
4th & 5th; (3) the region 7560-7730A for Oct. 9th & 10th; (4) the region 7000-7220A for Oct. 9th & 10th. The dashed lines are the best 
fit to each panel (see Table Eft. The full lines are the weighted mean of the individual fits. 




7600 7700 7B00 7900 8000 8100 8200 



Wavelength in Angstroms 

Figure 7. Phase-binned spectra of IP Peg from 1995 October 4th 
&; 5th. Phase ranges are marked, as are the Kl and Nal absorption 
features. 



used for the cross-correlation, which is presumably why they 
have no significant effect on the radial velocity study. 



5 A NEW EPHEMERIS 

Our spectroscopically derived eclipse timing-point and that 
of M89, both corrected to time of mid-white dwarf egress, 
were added to the photometric timings of Wood et al. (1989) 
and those of Wolf et al. (1993). All timings were converted 
to TDB. A linear ephemeris was then calculated to give 

JD(TDB) T ogrcss = 244 7965.891 44(5) + E * 0.158 206 17(3) 

which is for the time of mid-white dwarf egress. Fig. ^ shows 
the O — C plot for this ephemeris. We note in passing that 
the O — C residuals calculated from the linear ephemeris can 
be fitted by a sine wave with parameters close to those given 
in Wolf et al. (1993), which they interpreted as a third body 
in the system. 
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Figure 8. The O-C plot for the mid-white dwarf egress 
ephemeris from the high-speed photometry timings of Wood et 
al. 1989 (circles), Wolf et al. 1993 (crosses) and the spectroscopic 
points from M89 and this paper (filled triangles, also marked with 
M and B for clarity). 



6 THE APPARENT ECCENTRICITY 

6.1 Are the new data consistent with a circular 
orbit? 

M89 found that an elliptical orbit was a better fit to their 
radial velocity data than a circular one. Here, we test our 
spectra for evidence of eccentricity. 

We fitted each red-star radial-velocity data set with a 
circular orbit of the form 

V = 7 + K 2 sin(27r[0 i -0 O ]) 

and then an elliptical orbit of the form (e.g. Smart 1977, 
Chap. XIV; Petrie 1962) 

V — 7 + K2 [cos(i/ + lo) + e cos(oj)] 

recording the mean square deviations (S 2 ) of the data from 
the fitted model in each case. In the elliptical orbit v is the 
true anomaly of the red star, u> is the longitude of periastron 
of the red star from the centre of gravity of the system, and 
e is the eccentricity of the red-star orbit. Following Lucy and 
Sweeney (1971), we use the F-test 



F 



Nobs — /ell 
fell /crc 



era 

'-'crc 
°cll 



(1) 



to determine the significance of the elliptical fit over the 
circular one. iVobs is the number of observations, / c n is the 
number of degrees of freedom in the elliptical fit, / crc those 
in the circular fit, and S% IC ,Seii are the variances of the 
data from the circular and elliptical fits respectively, i.e. 
S 2 = T,(ydata — Ufit) 2 - This test may be derived from the 
relationship between the F- and x 2 - distributions (e.g. Bev- 
ington 1969, Ch. 10). For both the circular and elliptical 
fits the systemic velocity, 7, was held fixed at zero since our 
method of radial velocity extraction uses an orbital-average 
as a template; the orbital period was also held fixed. The 
number of free parameters in each fit is then / e u = 4 (K2 ,e,v 
and u>) and / cr c = 2 (K2 and (f>o). 

In statistical terms we are testing the hypothesis that 



the radial velocity curves are consistent with zero eccen- 
tricity. Table |^ shows for each radial velocity data set the 
probability (p) that this hypothesis is correct, and since p 
never falls below 30 percent, we must accept the hypothesis. 



6.2 Are the new data consistent with the 
previously measured eccentricity? 

Although our data are consistent with a circular orbit, it is 
still possible that, were they noisy enough, they could also 
be consistent with an eccentricity at the level measured by 
M89. So we determined the maximum level of eccentricity 
that could be present in our data. 

For this exercise we fitted simultaneously the veloci- 
ties from both TiO bands, but retained the separate data 
sets for each resolution. For each data set we constructed 
an e-\ 2 space by freezing e at set values but allowing all 
other parameters (except orbital period, always held fixed) 
to run free. Imposing the condition that xt~^ f° r the e~0 
fit, defines the error for each velocity point. From the 90% 
confidence interval for five free parameters (Lampton, Mar- 
gon & Bowyer 1976), both data sets yield an upper limit to 
any eccentricity in our data of around 0.05. 

Thus we conclude that our upper limit to the eccentric- 
ity is smaller than M89's detection (0.089±0.020). 



6.3 Is the difference in eccentricity due to a 
difference in technique? 

There are two important differences between our measure- 
ments and those of M89. They used the Nal doublet, whilst 
we used TiO; and they cross correlated against an M-star 
template, whilst we used a velocity-shifted mean of the data. 
So, we performed a cross-correlation study using the M3.5V 
template GL27B, over the Nal region used by M89. For test- 
ing the significance of the eccentricity we used / e n = 5 and 
/crc = 3, since 7 is an extra free parameter in each case. The 
results are shown in Table ^| and show that this study also 
returns zero eccentricity. 

The only remaining difference in technique is that M89 
used a slightly different formula for the radial velocity, which 
breaks down for e > 0.1. Although this approximation 
should not create a problem, for completeness we also tested 
their table 2 data, labelled M89 in our Table |^. Our code 
finds a similar eccentricity (e = 0.094), rejecting a circular 
orbit at greater than 99.99% confidence. 



6.4 Orbital eccentricity; a discussion 

The conclusion is thus inescapable, the apparent eccentric- 
ity of the orbit has changed between the observations of 
M89 and our own. IP Peg is not the only system where an 
apparent eccentricity has been observed in the radial ve- 
locity data. AM Her, U Gem, CH UMa and YY Dra all 
show significant eccentricities (Friend et al. 1990a, b). The 
explanation that has been put forward to explain these de- 
viations is irradiation of the secondary star, either by the 
white dwarf, bright spot or a combination of both. Under 
this assumption, Davey & Smith (1992) were able to map 
the distribution of the Nal doublet across the surface of the 
secondary star. As Davey & Smith point out, there is no 
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satisfactory theoretical explanation for the asymmetric dis- 
tributions this model predicts, nor is there any apparent 
correlation with other parameters which would indicate the 
presence of strong heating. As our spectra show the Nal dou- 
blet to be strongly contaminated by emission from the ac- 
cretion disc, we put forward an alternative hypothesis, that 
this contamination, even when it is much weaker, affects the 
radial velocity curves of CVs. 

IP Peg is a stunning example of the random nature 
of deriving an eccentricity from radial velocity data: it has 
been observed twice, over the same wavelength region, both 
times in quiescence (although the M89 data may have been 
on the decline from a small outburst), both use a similar 
procedure, both return the same value for the radial velocity 
of the secondary star and yet one shows an eccentricity while 
the other does not. The raison d'etre for M89's correction to 
the radial velocity was the apparent eccentricity of the orbit. 
Significant irradiation would introduce an eccentricity (see 
M89), so the the absence of such effects (see also Section 
|i|), means that we should not make any correction to our 
semi- amplitude . 



7 COMPONENT MASSES 



The mass function 



/(M) 



2tvG ' 



(2) 



where the symbols have their usual meaning, links the dy- 
namical motion of the secondary star to the mass of the 
unseen component. Using the values obtained in Section |§| 
we find 



/(M) = 0.596 ±0.031 M s 



(3) 



which acts as a lower limit to the mass (Mi) of the white 
dwarf star in IP Peg. Using this value, we can calculate 
masses of the individual components if the mass ratio, q = 
M2/M1 (where Mi is the mass of the white dwarf), and 
inclination, i, are known using 



Mi = (1 



. f(M) . 
sin 3 (z) ' 



M 2 = qM Y ■ 



(4) 



The mass ratio may be determined by combining our K2 
with the rotational broadening measurement of Catalan, 
Smith & Jones (2000). They find V ro t sin(i) =125±15 km 
s _1 where 



Vrot sin(i) 
~K 2 



= (! + «)/(<?) 



(•>) 



and f(q) is any of the approximations to the fractional 
Roche- lobe radius Rija. Here we use the approximation 
due to Eggleton (1983), which gives a mass ratio of q — 
0.322 ± 0.075. We derive the inclination using the eclipse 
width of A<f> = 0.0863 from Wood & Crawford (1986). To 
be consistent with this eclipse width, we must raise our lower 
limit on q to 0.285, since smaller mass ratios will not produce 
a long enough eclipse, even if i = 90 degrees. This sets the 
upper bound of the inclination to 90 degrees, and thus our 
measurement of the inclination is i = 85.9^2'g degrees. Since 
the error in V ro t sin(i) dominates the errors for q, i, Mi and 
M2, we can simply use the extreme values of i, with their 



corresponding values of q in equation Q to obtain our lim- 
its on the component masses. These are Mi = l.O5±°;J*M 
and M 2 = O.33tj};^M . 

Martin et al. (1987) find, from the reduction in cross- 
correlation errors, that the secondary star is probably closer 
to M4.5V than M3.5V. From the TiO line-depth ratio, 
Catalan et al (2000) find the spectral type to be M5±0.5V. 
Using the observationally derived mass-spectral type rela- 
tion of Kirkpatrick & McCarthy (1994-their equation (11)), 
which is supported by the theoretical models of Baraffe & 
Chabrier (1996), our value of M 2 = 0.35 ± 0.08M© implies 
a spectral type of M2.2±0.6V. The error in the calculated 
spectral type includes both the error in our mass estimate 
and the error in the relation itself. For a spectral type of 
M4V, this same relation returns a mass of O.I8M0. Although 
our mass estimate is now consistent with an M-dwarf, the 
secondary star in IP Peg is still over-massive for its derived 
spectral type. 

We can also calculate the radius for the secondary star, 
and for consistency do so using equation (3) of Smith & 
Dhillon (1998). This gives a radius of R 2 = 0.40 ± O.O3R . 
We can then place the star in the mass-radius-period re- 
lationships given in that paper. In all cases IP Peg is now 
consistent not only with the other CVs, but also with the 
current semi-empirical and theoretical relationships. 

Our value for q is low compared with previous determi- 
nations, and although the error bars overlap with those of 
Wood & Crawford (1986) who find 0.35 < q < 0.49 from 
measurements of the bright spot, it is probably inconsistent 
with that of Wolf et al (1993) who measure q = 0.6 (no er- 
ror given) from a similar technique. Both values rely on the 
center-of-light for the stream following a ballistic trajectory, 
which may well not be the case once it begins to interact 
with the Keplerian flow in the disc. Marsh (1988) also find 
q = 0.6, but his value relies on measuring the radial veloc- 
ity of the primary using the wings of emission lines from the 
disc. Long & Gilliland (1999) have shown how this technique 
is inaccurate for U Gem. 



8 THE SECONDARY STAR LIGHTCURVE 

Measuring the flux in the TiO bands should allow us to ex- 
tract the flux of the secondary star alone, assuming there 
is no contribution to the TiO from the other sources of the 
light in the system. Wade & Home (1988) constructed such 
a lightcurve for the short period dwarf nova Z Cha, by mea- 
suring the total flux in the TiO band with respect to a fitted 
continuum (i.e. a flux deficit). Here we undertake a similar 
study, but using a slightly more sophisticated method. 



8.1 The method 

If the temperature dependence of the equivalent width of the 
spectral features is small for the range of temperatures over 
the secondary star, the flux from the star at any given phase 
is some phase-averaged secondary-star spectrum multiplied 
by a scaling factor we refer to as k^. Thus, if we scale and 
then subtract the phase-averaged spectrum from the spec- 
trum at a given phase, the residual will be smooth across the 
TiO bands when we have chosen the right k^. Since both the 
phase- averaged spectrum, and the spectrum at the phase 
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Figure 9. The fc^-factor plotted against phase for the lower 
resolution data ( Oct. 4th & 5th, upper curve, displaced upwards 
by 0.25 for clarity) and for the higher resolution data (Oct. 9th & 
10th, lower curve, displaced downwards by 0.4). The region fitted 
was 7450-7750A. Shown over-plotted (full line) is the best-fitting 
'ellipsoidal modulation with eclipses' model. The data are plotted 
twice and k$ is normalised to be about one at phase zero. 
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Figure 10. As the previous figure but derived fitting the region 
7040-7340A in the higher resolution data. 

under consideration have contributions from the accretion 
disc, the residual is not simply the accretion disc flux at the 
phase in question. However, if the disc has a spectrum which 
is approximately a straight line at all phases, the residual 
will also be a straight line. 

To determine fc<p for a given spectrum we subtract from 
it varying portions of the mean spectrum, fit the residual 
with a straight line, and evaluate \ 2 f° r the fit. We then 
look for the value of k^ which minimises x 2 ■ To minimise 
the effect of the assumption of a straight line fit, we evaluate 
\ 2 over a relatively short fragment of spectrum. 

8.2 Results 

Fig. [] shows kj, plotted against orbital phase calculated 
from the region between 7450-7750A for both the high and 
low resolution data sets. The average (x 2 ) s (over all spec- 



tra) for the individual straight-line fits to the residuals were 
2.21T0.23 and 3.59T0.45 for the high and low resolution 
data, respectively (the error is the standard deviation of the 
individual fits around these means). Fig. |l(j shows k^ plot- 
ted against orbital phase for the case of a straight line fit to 
the region between 7040-7340A for the low resolution data 
set. 

The most prominent feature in Figs. ^ and [lo] is the 
dip at phase 0.5. The obvious interpretation for this is that 
it is an eclipse of the secondary star by the accretion disc. 
This is a rather remarkable result, since the accretion discs 
of quiescent dwarf novae are normally thought to be opti- 
cally thin. However, we are viewing the disc in IP Peg at an 
unusual angle (nearly edge on), and whilst it may be that 
the discs are optically thin when viewed from the pole, the 
large column length to the secondary star at phase 0.5 may 
be enough to obscure it. Alternatively, it may be that there 
are some optically thick regions in the disc, which cause the 
eclipse we see. 

At other phases the lightcurves are broadly consistent 
with the double-humped ellipsoidal modulation we might ex- 
pect. To see just how close they are to this expectation, and 
prove that the suggestion of an eclipse is at least plausible, 
we modelled the data using the code described in Ioannou 
et al (1999). The model is a standard ellipsoidal model with 
Roche geometry, but importantly for this work, includes the 
mutual eclipses between the disc and secondary star. Here 
the disc is assumed to be a cold (i.e. dark) obscurer. We used 
the parameters derived in Section |?] in addition to a limb 
darkening coefficient of 0.4, a gravity darkening exponent of 
0.08 and a secondary star (pole) temperature of 3375K. We 
performed a grid search in opening angle and radius to find 
the best fit to the data. 

In Figs. ^ and [l(] we plot the best fitting models. This 
shows that the model can account for the broad outline of 
the lightcurves, but there are problems in detail. This is 
reflected by the fact that the best xt obtained was 6.1 (for 
the data of Fig |l^), and that all the datasets required a 
disc radius equal to the Roche-lobe radius of the primary 
(the maximum we allowed) . For the lower resolution data of 
Fig. ^ the best fit opening angle was 8°, 4° for the higher 
resolution. For Fig. holthe opening angle was 3°. 

The models throw into sharp relief how the flux at (j> = 
0.75 is less than that at cj> — 0.25. Furthermore, Fig. [] shows 
how this asymmetry can vary from night to night. Similar 
night-to-night variability was seen in the Z Cha TiO flux 
deficits by Wade & Home (1988). The asymmetry, and its 
variability, could be explained by strong irradiation from 
the bright spot depleting the TiO on the side of the red 
star facing the bright spot, i.e. the side seen at <fi = 0.75. 
However, this would produce a marked eccentricity in the 
radial velocity curves deduced from the TiO regions of the 
spectrum. As shown in Section ^, no significant eccentricity 
is found from any of the IP Peg absorption features. 

The models suggest that the problem is with the data 
around phase 0.25. This means obscuration of the secondary 
star by the disc is unlikely, given the star's position well clear 
of the disc at this phase. Thus we suspect the problem lies 
with our simple contamination model (a straight line fit). 
The spectral variations of the disc with time are probably 
more complex. 
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9 CONCLUSIONS 

The main conclusions from this work are as follows. 

(1) There is some form of variability which means that radial 
velocity studies of the infrared Nal doublet in IP Peg can 
sometimes return an apparently elliptical orbit. We suspect 
this is contamination from disc emission. 

(2) We strongly recommend the use of TiO rather than Nal 
for future radial velocity studies. Not only does it avoid the 
contamination problem described above, but it gives a fac- 
tor two gain in signal-to-noise. 

(3) Our TiO study suggests the radial velocity semi- 
amplitude for the secondary star should be revised to K2 = 
331.3 ± 5.8km s _1 , and as a result the primary and sec- 
ondary star masses become Mi = 1.05lg J7M0 and M2 = 
0.33l°' J5M0 respectively. 

(4) Despite this downwards revision of the secondary star 
mass, it is still over-massive for its observed spectral type, 
but now agrees with the current semi-empirical and theoret- 
ical mass-radius-period relationships. 

(5) The accretion disc eclipses the secondary star in quies- 
cence, implying that, when viewed at high inclination, it is 
optically thick. 
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